Spin-orbit mediated anisotropic spin interaction in interacting electron systems 
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We investigate interactions between spins of strongly correlated electrons subject to the spin-orbit 
interaction. Our main finding is that of a novel, spin-orbit mediated anisotropic spin-spin coupling 
of the van der Waals type. Unlike the standard exchange, this interaction does not require the wave 
functions to overlap. We argue that this ferromagnetic interaction is important in the Wigner crystal 
state where the exchange processes are severely suppressed. We also comment on the anisotropy of 
the exchange between spins mediated by the spin-orbital coupling. 



PACS numbers: 71.70.Ej, 73.21.La, 71.70.Gm 

Introduction. Studies of exchange interaction between 
localized electrons constitutes one of the oldest topics in 
quantum mechanics. Strong current interest in the pos- 
sibility to control and manipulate spin states of quantum 
dots has placed this topic in the center of spintronics and 
quantum computation research. As is known from the 
papers of Dzyaloshinskii [l[ and Moriya [2| , in the pres- 
ence of the spin-orbital interaction (SOI) the exchange is 
anisotropic in spin space. 

Being a manifestation of quantum tunneling, the ex- 
change is exponentially sensitive to the distance between 
electrons Q. This smallness of the exchange leads to 
a large spin entropy of the Wigner crystal state, as 
compared to the Fermi liquid state, of diluted two- 
dimensional electron gas in semiconductor field-effect 
transistors Qj- The consequence of this, known as the 
Pomeranchuk effect, is spectacular: Wigner crystal phase 
is stabilized by a finite temperature. 

In this work we show that when subjected to the 
spin-orbit interaction, as appropriate for the structure- 
asymmetric heterostructures and surfaces [sl , interacting 
electrons acquire a novel non-exchange coupling between 
the spins. The mechanism of this coupling is very similar 
to that of the well-known van der Waals (vdW) interac- 
tion between neutral atoms. This anisotropic interaction 
is of the ferromagnetic Ising type. It lifts extensive spin 
degeneracy of the Wigner crystal and leads to the long- 
range ferromagnetic order. We also re- visit and clarify 
the role of spin-orbit interaction in lowering the sym- 
metry of the exchange coupling between spins. Particu- 
larly, we point out that the exchange Hamiltonian, de- 
spite its anisotropic appearance, retains spin-rotational 
invariance to the second order in the spin-orbital cou- 
pling. We argue that spin-rotational symmetry is broken 
only in the forth order in SOI coupling. 

Calculation of the vdW coupling. To illuminate the 
origin of the vdW coupling, we consider the toy prob- 
lem of two single-electron quantum dots described by the 
double well potential [1, Q , see Figure [U 
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FIG. 1: (color online) Two-dot potential ([T}. Blue (dark grey) 
arrows indicate electron's spins. 

where uj^^y are confinement frequencies along x/y direc- 
tions. The electrons, indexed by j = 1,2, are subject to 
SOI of the Rashba typefHI with coupling an 

Hso ^ ^ oiR Vj X ■ ^, (2) 

where Si are the Pauli matrices and z is normal to the 
plane of motion. Finally, electrons experience mutual 
Coulomb repulsion so that the total Hamiltonian reads 

^ 2 

H=Y.^^ + ^^^i^yj)\ + + ^^so. (3) 

At large separation between the two dots the exchange 
is exponentially suppressed and the electrons can be 
treated as distinguishable particles. One then expects 
that Coulomb-induced correlations in the orbital motion 
of the electrons in two dots translate, via the spin-orbit 
interaction, into correlation between their spins. Con- 
sider the distance between the dots, a, much greater 
than the typical spread of the electron wave functions, 
l/^TOWj;. In this limit the electrons are centered about 
different wells, and the potential can be approximated as 

V{ri,r2) « \mul{{xi - a/2f + {x2 + a/2f) 

+ \'n^^l{vl + vl)- (4) 
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At this stage it is crucial to perform a unitary transfor- 
mation Q which removes the hnear spin-orbit term 
from © 



U = exp[?mai?z • (ri x (Ti + r2 x <?2)]- 



(5) 



Owing to the non-commutativity of Pauh spin matrices, 
SOI can not be ehminated completely, resulting in higher 
order in the Rashba coupling contributions as given 
by iJ = UHW below 



(6) 



Here Lj is the angular momentum of the j electron, 

= xpy — ypx, and tilde denotes unitarily rotated opera- 
tors. The calculation is easiest when the confining energy 
is much greater than both the Coulomb energy e^/a and 
the spin-orbit energy scale y/mcoaji. In terms of the new 
(primed) coordinates r'l = ri — d/2 and r'2 = r2 + a/2 
centered about (a/2,0) and (—a/2,0), respectively, the 



interaction potential e^/|r'i — r'2 



is expanded in 



powers of 1/a keeping terms up to second order in the 
dimensionless relative distance (r'l — r'2) /a- The linear 
term, e'^{xi — X2)/a^, slightly renormalizes the equilib- 
rium distance between the electrons and can be dropped 
from further considerations. In terms of symmetric (S) 
and anti-symmetric (A) coordinates: 



XS/A 



e'i lb x'2 

~7T~ 



VS/A 



^/2 ' 



(7) 



the quadratic term e^{2{x'^ 



renormalizes the anti-symmetric frequency luax ~* + 
4e^/(ma'^) and uj\ — ui'^ — 2e^/{ma'^) while leaving the 



ljI and w|y = u)y. 



symmetric ones unmodified, ujsx'^ 

Quite similarly to the textbook calculation of the vdW 
force jlOj], the resulting Hamiltonian H = Hg+Hji+HsQ 
becomes that of harmonic oscillators 



Hs/A — 



Ps/A 

2m 



+ y {^Is/A^l/A + ^Is/Ayl/A) (8) 



perturbed by i?so — H^q + SHqq + 0{a'j^), where 



(2) 



(2) 
SO 



ma. 



■[{xsPyS - ysPxs) + S ^ A]{al + dl) 



^ [{xsPyA - yApxs) + S" ^ A] (ctJ - a| ) , (9) 



"-"so 



2^/2 



'■[Pysi^l - Ol) + PyAi^i + i^l)]- (10) 



It is evident from Eqn. (lOTTOl) that the leading correc- 
tions to the ground state energy is obtained either by the 
excitation of a single y-oscillator (through ^TU^ ) and by 



the simultaneous excitation of oscillators in both the x 
and y directions (through ([9])), 



i,j=S,A 



|(0|5i/('^|ly,)P , mH^^l\lx,ly,)\' 



Jiy 



^3V 



It is easy to see that the spin-dependent contributions 

from SHgQ cancel exactly while those originating from 
~ (2) 

HgQ do not, resulting in the novel spin interaction 

i?vdW = ^m'^aj^ a-ld-2(4'{uJSy,UJSx) + 4>{i^Ay,0JAx) 



(uJAy, l^Sx) - 4>[l-^Sy, t^Ax)^ , 



(11) 



where the function <f> is given by a simple expression 



<t){x,y) 



{x-yf 
xy{x + y)' 



In case of cylindrically symmetric dots. 



vdW 



4a6, 



a^a2- 



(12) 



(13) 



The physics of this novel interaction is straightforward: it 
comes from the interaction-induced correlation of the or- 
bital motion of the two particles, which, in turn, induces 
correlations between their spins via the spin-orbit cou- 
pling. The net Ising interaction would have been zero if 
not for the shift in frequency of the anti-symmetric mode 
due to the Coulomb interaction. Note that the coupling 
strength exhibits the same power-law decay with distance 
as the standard van der Waals interaction [iq|. 

From (fTTjl . it follows that in the extreme anisotropic 



limit of 



cxD, or equivalently, the one-dimensional 



(ID) limit, there is no coupling between spins. This result 
is understood by noting that ID version of SOI, given by 
Q!fl 'Ylij "'J-Pj ; can be gauged away to all orders in aji by a 
unitary transformation Uxb = ex^\imafi(x\G\ -|-X2cr|)]. 
Hence the absence of the spin-spin coupling in this limit. 
However, either by including magnetic field (Zeeman in- 
teraction, see below) in a direction different from cr^, or 
by increasing the dimensionality of the dots by reducing 
the anisotropy of the confining potential, the spin-orbital 
Hamiltonian acquires additional non-commuting spin op- 
erators. The presence of the mutually non-commuting 
spin operators (for example, cr^ and ct^ in (l2|)) makes it 
impossible to gauge the SOI completely, opening the pos- 
sibility of fluctuation-generated coupling between distant 
spins, as in equation (|13[) . 

Effect of the magnetic Held. For simplicity, we neglect 
orbital effects and concentrate on the Zeeman coupling, 
Hz = — J2j /'^^ where = g^B- Unitary transfor- 
mation ([5]) changes it to Hz~Az'maRa{af-'a2)/2+SHz- 
Here 



6Hz 



E 

i=i,2 



TOQ!/?A^(x'o-^ 



(14) 
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describes the couphng between the Zeeman and Rashba 
terms. In the basis (O it reduces to 



SHz = -mA^aR 



V2 



-mA^aR -j= . (15) 

For sufficiently strong magnetic field, 3> ^/mcoaR, 
Hso can be neglected in comparison with 6Hz- Calcu- 
lating second order correction to the ground state energy 
of the two dots, represented as before by Hs + Ha, and 
extracting the spin-dependent contribution, we obtain 



AE^ 



-AW 



(2- 



1"2 



(16) 



In the extreme anisotropic limit ujy 00 the dots become 
ID and we recover the result of Ref. JJ,. For the isotropic 
limit LOx — ^y, the coupling of spins acquires a rnagnetic 
dipolar structure identical to that found in Ref. 112 . 

Anisotropy of the exchange. Next, wc allow for the 
electron tunneling between the dots. The spin dynamics 
of the electrons is now described by the sum of exchange 
and the van der Waals interactions, H = H^x + H^dw- 
Here the exchange coupling, H^xi contains both isotropic 
and possible anisotropic interactions, while H^dw is given 
by (fTTj) and (|13|1 . In the absence of spin-orbit interac- 
tion, the total spin is conserved and the Hamiltonian is 
SU(2) invariant. As such, the only spin interaction al- 
lowed has the well known isotropic form i?Ex (?i • (T2. 
The anisotropy of the exchange is mediated by the spin- 
rotational symmetry breaking SOI 0. When the tun- 
neling is no longer spin-conserving, electron spins precess 
while exchanging their respective positions, giving rise to 
the anisotropic terms. As a result 13l. Il4l. llSj 



^^Ex = ^(bai-a2+Dd-(Ti x 0^2 +r(d-CTi)(d-CT2)) , (17) 

where d is the unit Dzyloshinskii-Moriya vector, of ampli- 
tude D, with odd dependence on the spin-orbit coupling 
aR. Coefficients b and F have even dependence on the 
spin-orbit coupling [3, [3] , while the exchange integral J, 
independent of aR, in this representation, sets the over- 
all energy scale. The direction of the DM vector can be 
understood as follows. As the D-term must be even un- 
der exchange operation P : 1 <-> 2, its amplitude must be 
odd with respect to inter-spin distance a = ri — r 2 = ax, 
hence d a = x. In addition, as z — > —z transformation 
in 12]) changes sign of aR, it must be that d ^ z as well. 
Thus, it must be that d = z x d — y. 

In the simplest approximation one neglects the "rem- 
nants" of SOI ([6]) altogether and writes the only possible 
exchange coupling H^J — -jCTi • (72 in terms of unitarily 
transformed spin operators aj . The meaning of this in- 
teraction is understood in the original basis by undoing 



the unitary transformation, iJg"^ = U'^ h'^^U . Using ^ 
and replacing ri, r2 by their respective average values, 
a/2x and —a/2x, one observes that spin 1 (2) is rotated 
about y axis by the angle = maRa in clockwise (coun- 
terclockwise) direction. As a result, one immediately ob- 
tains Eq. (jl7[) with parameters 



bo = cos 20, Do = sin 2d, Fq = 1 - cos 261, d = y. (18) 

As it originated from the S'J7(2)-invariant scalar product 
CTi • (72, the Hamiltonian (fTT)) with parameters (HH]) does 
not break spin-rotational SU(2) symmetry, despite its 
asymmetric appearance. Because of its "non-diagonal" 
nature, the D-term affects the eigenvalues only in ^ 
6^ order, and must always be considered together with 
the F-term. In the current situation psp , the two contri- 
butions compensate each other exactly. This important 
observation, made in Ref ill (see also [I^), was over- 
looked in several recent calculations of the DM term 

Uliilii. 

It is thus clear that the symmetry-breaking DM term 
must originate from so far omitted Hso fS]). To cap- 
ture it, we set up the exchange problem calculation along 
the lines of the standard Heitler-London (HL) approach. 
Despite its well-known shortcomings [20|, |2l|, 22], this 
approach offers conceptually simple way to estimate ex- 
change splitting and the structure of anisotropic spin 
coupling. Our basis set is formed by the antisymmetrized 
two-particle wave function \tp) — \tp) — P\ip), 

IV^) =^(l,2){ci| TT>+C2| Ti>+C3| iT>+C4| ii>} (19) 

is written in terms of unknown coefficients ci_4. Here 
(^(1,2) = f{xi - a/2)f{yi)f{x2 + a/2)/(y2) describes 
spatial wave function of distinguishable particles lo- 
calized near (a/2,0) and (—a/2,0), respectively, and 
f{x — xo) denotes the ground state wave function of one- 
dimensional harmonic oscillator centered around x — xo. 
As constructed, Lp{l,2) is the lowest energy eigenstate of 
two particles moving in the potential profile (|4]). 

The rest of the confining potential, Eq.([T]), together 
with the SOI (O, forms the perturbation 

^pert(l, 2) = ^ V{xj,y,) ~ Vin,f2) + Hso, (20) 

which is responsible for removing spin degeneracy of 
states contributing to (fT9|) . The eigenvalue problem 



{Ho + Vp,M^E\^), 



(21) 



where Ho is the sum of kinetic energy and confinement 
potential ([4]), is formulated as a 4 x 4 matrix problem 
by multiplying ([2T|) by the bra (siS2|</3(l, 2) from the left 
(here Sj=i,2 =T or J,) and integrating the result over the 
whole space. The obtained exchange Hamiltonian for the 
rotated spins a is of the form ([TT]) with 



2 



32maj^ 
9uJx<jJya ' 



(22) 
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while b = 1, r = to this order. The calculation sketched 
is valid in the large separation limit, a 3> l/y/muJ^^ and 
its most important feature is the scaling D ~ a'fj between 
the DM coupling and the spin-orbital one. This result is 
due to the fact that 0{aj^) term in ^ excites both x 



and y oscillators. Since the wave function contains 
only the ground states of the oscillators, the 0{aj^) term 
drops out and the first asymmetric correction originates 
in O(a^) terms of We checked that this crucial fea- 
ture is not an artifact of the HL approximation and is also 
obtained from a more reliable "median-plane" approach 
SEilli,!!!!, which we initiated. 

Noting that the DM term Dy • ai x a'2 affects the eigen- 
value of the two-spin problem only in order, we con- 
clude that exchange asymmetry due to the spin-orbit in- 
teraction may appear only in aj^ or higher order. This 
is because the effect of F-term in (|17l) on the eigenval- 
ues is of first order in F, and our calculation shows that 
F ~ 0(a|j,). Being proportional to J, see ((22)) . this con- 
tribution is also exponentially small. We then conclude 
that the leading source of spin anisotropy is provided by 
the vdW contribution (fTT|) and ([T3|) , which does not con- 
tain an exponential smallness of the exchange. 

Estimate of the vdW coupling. We now turn our at- 
tention to physical manifestations of the vdW spin cou- 
pling in the Wigner crystal. Neglecting the exchange 
interaction for the moment, we consider a two-electron 
problem within the frozen lattice approximation in which 
all other electrons are assumed fixed in their equilibrium 
lattice positions. The potential energy then is just that 
of four harmonic oscillators 2J] with frequencies w^^^ = 



a/(7=F 2)/{m'^aBa^) and lOu,v = \/ h ^ I) / {m^aBa^), in 
notations of Ref|2J. Here 7 « 5.52, [2J| as — K/{me^) 
is the Bohr radius, k is the dielectric constant and a 
is the lattice constant of the electron crystal, inversely 
proportional to the electron density n: a = (2/V3n)^/^. 
Repeating the steps that led to we obtain for the 
Wigner crystal problem 



vdW 



5vdW ^1^2 



(23) 



where B = Wt,) 4- (/)(cj^, ^.^n ) - (?!)(cj^,w„) - 

0(w^,a;„)]/8 — —3.75 • IQ^^y/rn^aBd^ ■ The spin-orbit 
mediated ferromagnetic coupling removes extensive spin 
degeneracy of the crystal, suppressing the Pomeranchuk 
effect physics [3]. Being of non-frustrated nature, it es- 
tablishes long-range magnetic order of Ising type with the 
ordering temperature of the order of the vdW constant 
<7vdw (|23|) . It should be compared with the much stud- 
ied Heisenberg exchange J„c = c(rs) exp[— 1.612 ^/r^, 
expressed in Rydbergs TZ — l/{2ma\). Here = 



1 / -y/ TTO^n is the dimensionless measure of the interac- 
tion strength, and the pre- factor c{rs) is a smooth func- 
tion of it 3] . We find that g^dw dominates the exchange 



for fs > r* « 20 in InAs, which has ap » 1.6 • lO^m/s 
[25! . For GaAs, with an « 300m/s more diluted 
situation is required, r* « 90. Given that multi-particle 
ring-exchange processes on the triangular lattice strongly 
frustrate any ordering tendencies due to the exchange [3|] , 
it appears that our estimate is just a lower bound on the 
critical density below which spin-orbit-induced ferromag- 
netic state should be expected. 
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